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^jrj' The classical trajectories of a particle governed by the PT-symmetric Hamiltonian 

^ . H = p^ + x^{ixY (e > 0) have been studied in depth. It is known that almost all 

\ trajectories that begin at a classical turning point oscillate periodically between this 

■<!:;;j- \ turning point and the corresponding "PT-symmetric turning point. It is also known 

that there are regions in e for which the periods of these orbits vary rapidly as 
' functions of e and that in these regions there are isolated values of e for which the 

classical trajectories exhibit spontaneously broken VT symmetry. The current paper 



pH , examines the corresponding quantum- mechanical systems. The eigenvalues of these 

^ \ quantum systems exhibit characteristic behaviors that are correlated with those of 

qh' the associated classical system. 
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O ■ I. INTRODUCTION 

od _ 

^ . This paper studies the eigenvalues of the non-Hermitian quantum systems defined by the 

PT-symmetric Hamiltonians 

- ^ . H = p^ + x^{ixy (e>0). (1) 

>< 

■ The Hamiltonians ([T]) can have many different spectra depending on the large- 1 a; | bound- 
ary conditions that are imposed on the solutions to the corresponding time-independent 
Schrodinger eigenvalue equation 

- ^"{x) + x'^{ixYi){x) = Ei){x). (2) 

The boundary conditions on ip{x) are imposed in Stokes' wedges in the complex-x plane. 
At the edges of the Stokes wedges both linearly independent solutions to ([2]) are oscillatory 
as — >■ oo. However, in the interior of the wedges one solution decays exponentially and 
the linearly independent solution grows exponentially. The eigenvalues E are determined by 
requiring that the solution ipi^x) decay exponentially in two nonadjacent wedges. Ordinarily, 
the eigenvalues are complex. However, if the two wedges are PT-symmetric reflections of one 
another, then there is a possibility that all the eigenvalues will be real. (The VT reflection 
of the complex number x is the number —x*.) 
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If e is not a rational number, there are infinitely many Stokes' wedges, which we label by 
the integer K. The center of the i^'th Stokes wedge lies at the angle 

_ (4ir + 2)7r _ vr 

['center ~ ^ _|_ ^ 2 ' ^ ' 

The upper edge of the wedge lies at the angle 

_ {AK + 3)7r _ vr 

Cupper ~ 4 _|_ g 2 



and the lower edge lies at the angle 



(4ir+l)7r TT 

f^lower — ^ _^ ^ 2' 

The opening angles of all wedges are the same: 

27r 

opening angle = ^ ^ . (6) 

The PT reflection of the wedge K is the wedge — i^' — 1. 

The "PT-symmetric eigenvalue problem posed in the wedges K = and K = —1 has been 
studied in depth in jl|-[3[. The spectrum for this problem is real when e > 0, as illustrated 
in Fig. [H For these boundary conditions the quantum theory specified by the Hamiltonian 
exhibits unitary time evolution [3]. 

The classical system corresponding to this quantum system has also been examined care- 
fully 0. Without loss of generality, one can take the classical energy to be unity. One 
can then graph the classical trajectories by solving numerically the system of Hamilton's 
differential equations 

dH . dH \ f ■ \f /n\ 

~dp ^ ^ " ~~dx' " ~ e)x{ix) (7) 

for an initial condition x(0); p{0) is determined from the equation H = E = 1. Since a;(0) 
and p{0) are not necessarily real and the differential equations ([7]) are complex, the classical 
trajectories are curves in the complex-x plane. For e > nearly all trajectories are closed 
curves. (When e is a positive integer, trajectories originating at some turning points can 
run off to infinity, but these are special isolated cases. When e < 0, all trajectories are open 
curves.) If e is noninteger, there is a branch cut in the complex-x plane, and to be consistent 
with VT symmetry we take this cut to run from to oo along the positive-imaginary axis. 
A closed trajectory may cross this branch cut and visit many sheets of the Riemann surface 
before returning to its starting point. 

As e increases from 0, the harmonic-oscillator turning points at a; = 1 (and at x = — 1) 
rotate downward and clockwise (anticlockwise) into the complex-x plane. These turning 
points are solutions to the equation 1 + (ix)^^'^ = 0. When e > is irrational, this equation 
has infinitely many solutions; all solutions lie on the unit circle and have the form 

—— g^^turning point 

where 

_ {2K+l)n _ n 

^turning point 2 + e 2 ' 
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FIG. 1: Eigenvalues of the Hamiltonian H m These eigenvalues are found by solving the 
Schrodinger equation ([2]) with boundary conditions imposed in the K = and K = —1 wedges. 
The eigenvalues are real, positive, and discrete when e > 0. When e = the spectrum consists of 
the standard harmonic-oscillator eigenvalues £'n = 2n + 1 (n = 0, 1, 2, . . .). As e increases, the 
eigenvalues grow and become increasingly separated. 

These turning points occur in PT-symmetric pairs (pairs that are symmetric when reflected 
through the imaginary axis) corresponding to the K values {K = —1, K = 0), {K = 
-2, K = 1), {K = -3, K = 2), {K = -A, K = 3), and so on. When e is rational, there 
are only a finite number of turning points in the complex-x Riemann surface. For example, 
when e = y , the Riemann surface consists of five sheets and there are eleven pairs of turning 
points. 

The period T of a classical orbit depends on the specific pairs of turning points that 
are enclosed by the orbit and on the number of times that the orbit encircles each pair. 
As explained in Refs. [1, [6|, any orbit can be deformed to a simpler orbit of exactly the 
same period. This simpler orbit connects two turning points and oscillates between them 
rather than encircling them. For the elementary case of orbits that oscillate between the 
K = —1, K = pair of turning points, the period of the closed orbit is a smoothly decreasing 



4 



Ti.^.2v-.mm^cos(^]. (10) 



function of e: 



r[(4 + e)/(4 + 2e)] V4 + 2e 

To derive ffTOl) we evaluate the contour integral §dx/p along a closed trajectory in the 
complex-x plane. This trajectory encircles the square-root branch cut that joins the pair of 
turning points. We deform the contour into a pair of rays that run from one turning point to 
the origin and then from the origin to the other turning point. The integral along each ray 
is a beta function, which is a ratio of gamma functions. Equation flTU]) holds for all e > 0. 

Complex trajectories that begin at turning points other than the K = Q and K = —1 
turning points are more interesting. For these trajectories it was shown in Refs. [Sj and 
that there are three qualitatively distinct regions of e. For the > turning point and the 
corresponding —K — 1 turning point, region I extends from e = to e = l/iiT. As e increases 
from in region I, the period of the classical trajectory is a smooth monotone decreasing 
function of e like that in (1101) . In region II, where e ranges from 1/K xvp to 4K, the period 
T(e) of the classical trajectories is a noisy function of e and the classical orbits in this region 
are acutely sensitive to the value of e. Small variations in e can cause huge changes in the 
topology and in the periods of the closed orbits. Depending on the value of e, some orbits 
have short periods and others have long and possibly even arbitrarily long periods. In this 
region there are isolated values of e for which the orbits are not symmetric when reflected 
through the imaginary axis, and thus these orbits are not VT symmetric; such orbits are 
said to have spontaneously broken classical VT symmetry. In region III, e ranges from AK 
up to oo and the period of the orbits is again a smooth monotone decreasing function of e. 

The abrupt changes in the topology and the periods of the orbits for e in region II are 
associated with the appearance of orbits having spontaneously broken VT symmetry. In 
region II there are short patches where the period is a relatively small and slowly varying 
function of e. These patches are bounded by special values of e for which the period of 
the orbit suddenly becomes extremely long. Numerical studies of the orbits connecting the 
i^th pair of turning points have shown that there are only a finite number of these special 
values of e and that these values of e are rational js], Some special values of e at which 
spontaneously broken PT-symmetric orbits occur are indicated in Figs. [2] and |3] by short 
vertical lines below the horizontal axis. 

Refs. [sl, IgI explain how very long-period orbits arise. In order for a classical particle 
to travel a great distance in the complex plane, its path must slip through the forest of 
turning points. When the particle comes under the influence of a turning point, it usually 
executes a large number of nested U-turns and eventually returns back to its starting point. 
However, for some values of e the complex trajectory may evade many turning points before 
it eventually encounters a turning point that takes control of the particle and flings it back 
to its starting point. Refs. js], [oj conjectured that there may be special values of e for which 
the classical path manages to avoid and sneak past all turning points. Such a path would 
have an infinitely long period. We still do not know if such infinite-period orbits exist. 

Ref. jsj did not explain why the period of a closed trajectory is such a wildly fluctuating 
function of e. However, in Ref. j6[ it was shown that for special rational values of e the 
trajectory bumps directly into a turning point that is located at a point that is the complex 
conjugate of the point from which the trajectory was launched. This turning point reflects 
the trajectory back to its starting point and prevents the trajectory from being VT symmet- 
ric. Trajectories for values of e near these special rational values have extremely different 
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FIG. 2: Period of a classical trajectory beginning at the K = \ turning point in the complex-x 
plane. The period is plotted as a function of e. The period decreases smoothly for < e < 1 (region 
I). However, when 1 < e < 4 (region II), the period is a rapidly varying and noisy function of e. For 
e > 4 (region III) the period is once again a smoothly decaying function of e. Region II contains 
short subintervals where the period is a small and smoothly varying function of e. At the edges of 
these subintervals the period abruptly becomes extremely long. Detailed numerical analysis shows 
that the edges of the subintervals lie at special rational values of e that have the form p/g, where 
p is a multiple of 4 and q is odd. Some of these rational values of e are indicated by vertical line 
segments that cross the horizontal axis. At these rational values the classical trajectory does not 
reach the K = —2 turning point and the VT symmetry of the classical orbit is spontaneously 
broken. The periods of the classical orbits are extremely long near these values of e in region II, 
so in this graph the period in regions I and III appears relatively small. 



topologies and thus have periods that tend to be relatively long. This explains the noisy 
plots in Figs. |2] and [31 

We are not certain whether for each turning point there are a finite or an infinite number 
of special rational values of e for which the classical orbit has a broken VT symmetry. The 
data used to produce Figs. |2] and |3] is far from exhaustive and the bars underneath the 
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FIG. 3: Period of a classical trajectory joining (except when VT symmetry is broken) the K = 2 
pair of turning points. The period is plotted as a function of e. As in the K = 1 case shown in 
Fig. [21 there are three regions. When < e < ^ (region I), the period is a smooth decreasing 
function of e; when ^ < e < 8 (region II), the period is a rapidly varying and choppy function of e; 
when 8 < e (region III), the period is again a smooth and decreasing function of e. 



horizontal axes are only the known examples of broken VT symmetry. There are more such 
special rational values of e = |. 

Having reviewed what is known about the classical situation, we come to the principal 
subject of this paper, which is to study what happens at the quantum-mechanical level (that 
is, how the eigenvalues behave) as e increases through regions I, II, and III. In Sec. [ITl we 
compare the classical and quantum theories for two cases, K = 1 (and K = —2) and K = 2 
(and K = —3). Our numerical studies show that the eigenvalues exhibit characteristic 
changes in behavior that correspond closely with the changes in behavior of the classical 
trajectories. In Sec. IIIII we make some concluding remarks. 
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FIG. 4: Eigenvalues associated with the K = I and K = —2 pair of turning points plotted as 
functions of e. (The periods of the classical orbits that begin at the K = 1 turning point are shown 
in Fig. [2] for regions I - III.) When e = 0, the spectrum is that of the harmonic oscillator. As 
soon as e increases from 0, the eigenvalues begin to become complex, starting with the high-energy 
part of the spectrum. When e = 1, there are no more real eigenvalues. Real eigenvalues begin to 
reappear at e = 3, and at the isolated point e = 4 all eigenvalues are real. The spectrum is once 
again entirely real when e is greater than about 5.5. 



II. BEHAVIOR OF EIGENVALUES AS e PASSES THROUGH REGIONS I III 
A. Eigenvalues corresponding with K = 1 and K = —2 turning points 

In Refs. jsl and ^ the behavior of the classical trajectories was studied and the period 
of the classical trajectories was found to have a remarkably complicated dependence on e 
(see Fig. [2]). In Fig. H] we plot the eigenvalues for the corresponding eigenvalue problem 
for < e < 7. (By the corresponding eigenvalue problem, we mean the eigenvalues for 
the eigenvalue problem associated with the K = 1 and K = —2 Stokes wedges.) These 
eigenvalues were calculated by numerically integrating the Schrodinger equation from large 
|x| into the origin along the centers of the two Stokes wedges and then applying the appro- 
priate matching condition on iIj{x) and ip^x) at x = 0. At the origin il){x) is continuous, 
but the matching condition on iIj'{x) incorporates the angle of rotation between the two 
wedges. These matching conditions are reminiscent of the techniques used by M. Znojil in 
his discussion of "quantum toboggans" 0|. 

The eigenvalues in Fig. H] exhibit the following features: At e = the spectrum is that 
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of the harmonic oscillator, En = 2n + 1 (ra = 0, 1, 2, . . .), and the entire spectrum is real, 
positive, and discrete. As soon as e increases from 0, the eigenvalues begin to coalesce into 
complex conjugate pairs; the process of combining into complex-conjugate pairs starts at 
the high-energy part of the spectrum. (This disappearance of real eigenvalues qualitatively 
resembles that shown in Fig. [T]as e decreases below 0.) When e reaches the value 0.5, only 
one real eigenvalue (the ground-state eigenvalue) remains. This eigenvalue becomes infinite 
at e = 1. Between e = 1 and e = 3 there are no real eigenvalues at all. The real ground-state 
eigenvalue reappears when e increases past 3, and as e approaches 4, more and more real 
eigenvalues appear. At the isolated point e = 4 the spectrum is again entirely real and agrees 
exactly with the spectrum in Fig. [T] at e = 4. At the isolated point e = 5 the spectrum is 
once again entirely real, but just above and below this value of e the spectrum becomes 
complex again. Finally, as e reaches a value near e = 5.5, the spectrum is again entirely 
real, positive, and discrete and it remains so for all larger values of e. 

Note that some of the qualitative changes in Figs. [2] and H] occur at the same value of e. 
For example, the noisy fluctuations in Fig. |2] begin at e = 1, just as the real ground-state 
eigenvalue blows up, and the fluctuations stop temporarily at e = 3, just as the real ground- 
state eigenvalue reappears in Fig. HI The noisy classical fluctuations cease completely at 
e = 4, just as the real eigenvalues begin to reappear in Fig. HI 

The most important qualitative changes in Fig. H] appear to be caused by the classical 
turning points entering and leaving the Stokes wedges inside of which the eigenvalue problem 
is specified (see Fig. [5]). For the case K = 1 the upper and lower edges of the wedge are 
given by 

6'upper = ^ _^ ^ — 2 ^^'^ Slower = ^ _^ ^ — (H) 

At e = the K = 1 classical turning point lies in the center of the wedge. However, as e 
increases, the turning point rotates faster than the wedge rotates, and at e = 1 this turning 
point leaves the wedge, as shown in Fig. [5](a). There are no turning points in the wedge 
between e = 1 and e = 3, which corresponds with the total nonexistence of real eigenvalues. 
At e = 3 a turning point re-enters the wedge, as shown in Fig. Mjo), but this is the K = 2 
turning point and not the K = 1 turning point. The re-entry of a classical turning point is 
strongly correlated with the reappearance of real eigenvalues. 

When the turning point leaves the wedge, WKB quantization is no longer possible. Thus, 
the high-energy portion of the spectrum is no longer real. In fact, we believe that there are 
no eigenvalues at all, real or complex, when 1 < e < 3. 



B. Eigenvalues corresponding with K = 2 and K = —3 turning points 

For a classical trajectory beginning at the K = 2 turning point, the period as a function 
of e again exhibits three types of behaviors. The period decreases smoothly for < e < ^ 
(region I). When | < e < 8 (region II), the period becomes a rapidly varying and noisy 
function of e. When e > 8 (region III) the period is once again a smoothly decaying function 
of e. These behaviors are shown in Fig. [31 As in Fig. |2l the period in regions I and III is 
very small compared with the period in region II. The classical trajectory that begins at the 
K = 2 turning point terminates at the K = —3 turning point except when VT symmetry 
is spontaneously broken. Broken-symmetry orbits occur at isolated points in region II. 

Figure [6] displays the eigenvalues of the corresponding quantum problem. Observe that 
the eigenvalues are all real at e = (the harmonic oscillator spectrum), but they they begin 
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(a) 



(b) 





FIG. 5: Turning points entering and leaving the K = 1 wedge as e increases. As e increases past 
1, the K = 1 turning point, whose angle is indicated by a dotted line, leaves the wedge, as shown 
by the arrow [panel (a)]. Then, as e increases past 3, a new turning point enters the wedge [panel 



to become complex as e increases from 0. There are no real eigenvalues when e reaches the 
value 1/2, just at the onset of noisy fluctuations in Fig. |3l There are several isolated values 
of e at 2, 4, 5, 7, 8, and 9, where the spectrum is entirely real. There are also two regions, 
1/2 < e < 3/2 and 5 < e < 7, where there are no real eigenvalues at all. Between 2 and 3 
there appears to be a patch where the eigenvalues are all (or almost all) real. (This patch 
corresponds with a temporary disappearance of noise in Fig. [3]). Real eigenvalues begin to 
reappear at e = 9, and the spectrum appears to be completely real beyond e = 21/2. 

The regions in which there are no real eigenvalues are exactly correlated with turning 
points leaving the Stokes wedges. At e = 1/2 the K = 2 classical turning point leaves the 
wedge at the angle 6 = 37r/2. Then, at e = 3/2 the K = 3 classical turning point enters 
the wedge at the angle 37r/2. Again, at e = 5 the K = 3 turning point leaves the wedge at 
the angle n/2 and at e = 7 the K = 4 classical turning point enters the wedge at the angle 



Figures [2] and [3] for K = 1 and Figs. H] and [6] for K = 2 illustrate a general pattern that 
appears to hold for all K. For classical orbits that oscillate between the i^th pair of turning 
points, there are always three regions, region I for which < e < region II for which 

< e < 4:K, and region III for which AK < e. When e = 0, the turning points lie on the 
real axis in the centers of the Stokes wedges. As e increases, the turning points rotate faster 
than the wedges rotate and they rotate out of the wedges at e = At this point there is 
a patch where there are no real eigenvalues. As e continues to increase, the K + 1 turning 
point rotates into the Stokes wedge and the eigenvalues begin to become real. Eventually, 
this turning point leaves the wedge and there are again no real eigenvalues. In total, there 
are K regions where there are no real eigenvalues. Eventually, the 2K turning point enters 
the wedge at the end of region II, and the spectrum begins to become entirely real and 
remains entirely real in region III. 



(b)]. 



9 = n/2. 



III. CONCLUDING REMARKS 
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FIG. 6: The eigenvalues associated with the K = 2 and K = —3 pair of turning points plotted as a 
function of e. (The periods of the classical orbits that begin at the K = \ turning point are shown 
in Fig. [3] for regions I - III. When e = 0, the spectrum is that for the harmonic oscillator. As soon 
as e increases from 0, the eigenvalues begin to become complex, starting with the high-energy part 
of the spectrum, and when e = 1/2, there are no more real eigenvalues. Real eigenvalues begin 
to reappear at e = 3/2, and there is a region of e near 5/2 where the spectrum is entirely real. 
There are no real eigenvalues between e = 5 and e = 7. At the isolated values e = 2, 4, 8, and 9 all 
eigenvalues are real. The spectrum is once again entirely real when e is greater than about 10.5. 

There are many isolated points at which the eigenvalues are all real. For example, when 
K = 2, real eigenvalues occur at 2, 4, 8, and 9. Some of these points can be analyzed in 
terms of VT reflection. Specifically: (i) From 1^ we can see that at e = 2 the K = 2 wedge 
is centered at the angle Tvr/G. The VT reflection of this wedge lies at the angle — 7r/6. Thus, 
the K = 2 spectrum agrees exactly with the K = spectrum in Fig. [T] for this value of 
e. (ii) At e = 4 the K = 2 wedge is centered at 37r/4. The VT reflection of this wedge 
is centered at 7r/4, and under complex conjugation this angle becomes — 7r/4. Thus, the 
K = 2 spectrum is identical to the K = spectrum for this value of e. This configuration 
of wedges is shown in Fig. [71 (iii) At e = 8 the K = 2 wedge is centered at vr/S. Under 
complex conjugation, this wedge is now centered at the angle — 7r/3, which corresponds with 
the wedge for the K = case. Thus, once again, the spectra for these two problems are 
identical at this value of e. (Note that for these three situations e is integer valued, and 
hence there is no cut in the complex-x plane.) 

While quantum mechanics and classical mechanics are of course different theories with 
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FIG. 7: The Stokes wedges for e = 4 for the cases K = and K = 2. Because the wedges are 
complex conjugate pairs the spectra are identical for this value of e. 

different regimes of applicability, our investigation has shown that there are nonetheless 
some correlations between the two. Regular classical behavior seems to correspond to a 
completely real quantum spectrum, while chaotic classical behavior in the complex plane is 
associated with either a partially complex spectrum or a complete absence of eigenvalues. In 
the systems we have studied, the onset of rapid and irregular variation in the periods of the 
classical trajectories correlates precisely with the disappearance of all quantum eigenvalues 
as the ground-state energy goes to infinity, while the end of the noisy region is associated 
with the gradual reappearance of a completely real spectrum. Within the region of rapid 
variation some correlations can be observed, but they are less clear-cut and merit further 
analysis. 
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